Cascades, thermalization and eddy viscosity in helical Galerkin truncated Euler flows 
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The dynamics of the truncated Euler equations with helical initial conditions are studied. Tran- 
sient energy and helicity cascades leading to Kraichnan helical absolute equilibrium at small scales 
are obtained for the first time. The results of [Cichowlas et al. Phys. Rev. Lett. 95, 264502 
(2005)] are extended to helical flows. Similarities between the turbulent transient evolution of the 
ideal (time-reversible) system and viscous helical flows are found. The observed differences in the 
behavior of truncated Euler and (constant viscosity) Navier-Stokes are qualitatively understood us- 
ing the concept of eddy viscosity. The large scales of truncated Euler equations are then shown 
to follow quantitatively an effective Navier-Stokes dynamics based on a variable (scale dependent) 
eddy viscosity. 



The role played by helicity in turbulent flows is not 
completely understood. Helicity is relevant in many at- 
mospheric processes, such as rotating convective (super- 
cell) thunderstorms, the predictability of which may be 
enhanced because of its presence [1]. However helicity, 
which is a conserved quantity in the three dimensional 
Euler equation, plays no role in the original theory of tur- 
bulence of Kolmogorov. Later studies of absolute equi- 
librium ensembles for truncated helical Euler flows by 
Kraichnan [2] gave support to an scenario where in helical 
turbulent flows both the energy and the helicity cascade 
towards small scales [3], a phenomena recently verified 
in numerical simulations [4-6]. The thermalization dy- 
namics of the non-helical spectrally truncated Euler flows 
were studied in [7]. However, Kraichnan helical equilib- 
rium solutions were never directly observed in simula- 
tions. Note that the Galerkin truncated non- helical Euler 
dynamics was recently found to emerge as the asymptotic 
limit of high order hypcrviscous hydrodynamics and that 
bottlenecks observed in viscous turbulence may be inter- 
preted as an incomplete thermalization [8]. 

In this letter we study truncated helical Euler flows, 
and consider the transient turbulent behavior as well as 
the late time equilibrium of the system. Here is a short 
summary of our main results. The relaxation toward 
a Kraichnan helical absolute equilibrium [2] is observed 
for the first time. Transient mixed energy and helicity 
cascades are found to take place while more and more 
modes gather into the Kraichnan time-dependent statis- 
tical equilibrium. It was shown in [7] that, due to the ef- 
fect of thermalized small-scales, the spectrally truncated 
Euler equation has long-lasting transients behaving simi- 
larly to the dissipativc Navier-Stokes equation. These re- 
sults, obtained for non-helical flows, are extended to the 
helical case. The concept of eddy viscosity, as previously 
developed in [7] and [9], is used to qualitatively explain 
differences observed between truncated Euler and high- 
Reynolds number (fixed viscosity) Navier-Stokes. Fi- 



nally, the truncated Euler large scale modes are shown to 
quantitatively follow an effective Navier-Stokes dynam- 
ics based on a (time and wavenumber dependent) eddy 
viscosity that does not depend explicitly on the helicity 
content in the flow. 

Performing spherical Galerkin truncation at wave- 
number k max on the incompressible (V • u = 0) and spa- 
tially periodic Euler equation 9 t u+(u- V)u = — Vp yields 
the following finite system of ordinary differential equa- 
tions for the Fourier transform of the velocity u(k) (k is 
a 3 D vector of relative integers satisfying |k| < fc max ): 



d t u a (k,t) = -^■p Q(37 (k)^u /3 (p,t)w 7 (k- p,t), (1) 



where V~ a f3-y 



kfjP a 



■ + k~ ( P a [3 withP Q/3 = S a/ 3-k a kp/k 2 . 



This time-reversible system exactly conserves the en- 
ergy E = J2k E ( k ^) and helicity H = Efc#(M), 
where the energy and helicity spectra E(k, t) and H (k, t) 
are defined by averaging respectively i|u(k',i)| 2 and 
u(k', t) ■ uj(— k', t) (u; = V x u is the vorticity) on spheri- 
cal shells of width Ak = 1. It is trivial to show from the 
definition of vorticity that \H(k,t)\ < 2kE(k,t). 

We will use as initial condition uo the sum of two ABC 
(Arnold, Beltrami and Childress) flows in the modes k = 
3 and k = 4, 



u (x,y, z) 



i (3) 

4 ABC 



(x,y,z) 



u 



(4) 
ABC 



(x,y,z) (2) 



where the basic ABC flow is a maximal helicity station- 
ary solution of Euler equations in which the vorticity is 
parallel to the velocity, explicitly given by 



UabcO^^ 2 ) = p {[ B cos(ky) 



Csin(£;,z)] x- 

[A sin(fcx) + C cos(fcz)] y + 
[A cos(kx) + B s'm(ky)] z} . 



(3) 



The parameters will be set to A = 0.9, B = 1, C = 1.1 
and uo = {A 2 + B 2 + C 2 )- 1 / 2 (l/3 4 + 1/A 4 )- 1 / 2 . With 
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this choice of normalization the initial energy is E = 0.5 
and hclicity H = 3 x 4 x (3 3 + 4 3 )/(3 4 + 4 4 ) = 3.24. 

Numerical solutions of equation (1) are efficiently pro- 
duced using a pseudo-spectral general-periodic code [10] 
with 512 3 Fourier modes that is dealiased using the 2/3 
rule [11] by spherical Galerkin truncation at k max = 170. 
The equations are evolved in time using a second order 
Rungc-Kutta method, and the code is fully parallelized 
with the message passing interface (MPI) library. The 
numerical method used is non-dispersive and conserves 
energy and helicity with high accuracy. 

Fig.l shows the time-evolution of the energy and he- 
licity spectra that evolve from (2) compensated by fc 5 / 3 . 
The plots clearly display a progressive thermalization 



The final values of a and f3 (when total thermalization 
is obtained) corresponding to the initial energy and he- 
licity are a = 4.12 x 10 7 and /3 = 7695. Therefore the 
dimensionless number /3 2 k 2 /a 2 is at most of the order 
10~ 4 and equations (4) thus lead to almost pure power 
laws for the energy and hclicity spectra, as is manifest 
in Figl.d. Fig. 1 thus shows for the first time a time 
evolving helical quasi-equilibrium. 

In order to analyze the run in the spirit of Cichowlas 
et al. [7] we define fcth(t) as the wavenumber where the 
thermalized power-law zone starts. We define the ther- 
malized energy and helicity as 



E t h(t)= #(M); #th(t)= J2 H ^ t ) 



(5) 
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FIG. 1: Compensated energy (•••) and helicity spectra (xx x) 
with the predictions (4) in solid lines and (6) in dotted lines, 
a) t = 4.8. b) t = 7. c)t= 10. d) t = 19.8. 



where E(k,t) and H(t,k) are the energy and helicity 
spectra. 

The temporal evolutions of i^hj-ffth and fcth(t) are 
shown in Fig. 2. 
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similar to that obtained in Cichowlas et al.[7] but with 
the non zero helicity cascading to the right. 

The truncated Euler equation dynamics is expected 
to reach at large times an absolute equilibrium that is 
a statistically stationary gaussian exact solution of the 
associated Liouville equation [12]. When the flow has 
a non vanishing helicity, the absolute equilibria of the 
kinetic energy and helicity predicted by Kraichnan [2] 
are 

E(k) ^ 4 " • H(k) ki(3 8n U) 
1 >~ a l-/3 2 P/a 2 ' [ >~ a 2 l-f3 2 k 2 /a 2 1 ' 

where a > and /3fc max < a to ensure integrability. The 
values of a and (5 are uniquely determined by the total 
amount of energy and hclicity (verifying \H\ < 2k max E) 
contained in the wavenumber range [1, A: max ] [2]. 



FIG. 2: a) Temporal evolution of E t h ( — ) ,H t h (•-•) and fcth(t) 
(• • • ) normalized by their respective initial values. E to t =0.5, 
Htot = 3.24 and fc max = 170. b) Left vertical axis: tem- 
poral evolution of eth = d "' (**★) and Navier-Stokes en- 
ergy dissipation e = 2i^o X^=i x k 2 E(k) (•••). Right ver- 
tical axis: r]th = d "^ th (* * *) and NS helicity dissipation 
V = vo Efe™r k 2 H(k) (o o o). 

The values of a(t) and /3(t) during thermalization can 
then be obtained from E t h(t) , H t h{t) and fcth(t) by in- 
verting the system of equations (5) using ^k 2 lax <C 1. 

The Kraichnan prediction (4) for the high-fc part of 
the spectra are shown (in solid lines) in Fig. 1. The plot 
shows an excellent agreement with the prediction. 

The low-A: part of the compensated spectrum presents 
a flat zone that amounts to A: -5 / 3 scaling for both the 



3 



energy and helicity spectra. This /c _5//3 behavior was 
predicted by Brissaud et al. [3] in viscous fluids when 
there are simultaneous energy and helicity cascades. The 
energy and helicity fluxes, e and r\ respectively, determine 
the prefactor in the inertial range of the spectra: 

E(k) ~ e 2 / 3 fc- 5 / 3 , H(k) ~ rye- 1 / 3 fc- 5 / 3 , (6) 

Helical flows have been also studied in high Reynolds 
number numerical simulations of the Navier-Stokes (NS) 
equation. Simultaneous energy and helicity cascades 
leading to the scaling (6) have been confirmed when the 
system is forced at large scales [4-6] . 

The energy and helicity fluxes e and 77 at intermediate 
scales in our truncated Euler simulation can be estimated 
using the time derivative of the thermalized energy and 
helicity: eth = ^jf^ and 77th = , whose temporal 
evolutions are shown in Fig. 2. The predictions (6) for 
the low-fc part of the spectra are shown (in dotted lines) 
in Fig. 1. The plot shows a good agreement with the 
data. Note that Fig. l.a corresponds to t = 4.8, that is 
just after the time when both the maximum energy and 
helicity fluxes (to be interpreted below as "dissipation" 
rates of the non-thermalized components of the energy 
and the helicity) are achieved, see Fig. 2. In this way E t h 
and H t h determine the thermalized part of the spectra 
while their time derivative determines an inertia! range. 

We now compare the dynamics of the truncated Euler 
equation with that of the unforced high-Reynolds num- 
ber NS equation (i.e. Eq.(l) with — v k 2 u a (k, t) added in 
the r.h.s.) using the initial condition (2). The viscosity is 
set to vq = 5 x 10~ 4 , the smallest value compatible with 
accurate computations using fc max = 170. A behavior 
qualitatively similar to that of the truncated Euler equa- 
tion is obtained (see Fig. 2b). However, the maxima of 
the energy and helicity fluxes (or dissipation rates) occur 
later, and with smaller values. 

We refered above to "dissipation" in the context of 
the ideal (time-reversible) flow. A proper definition of 
dissipation in the truncated Euler flow is now in or- 
der. Thermalized modes in truncated Euler are known 
to provide an eddy viscosity v c ddy to the modes with 
wavenumbcrs below the transition wavenumber [7]. It 
was shown in [9] that Monte-Carlo determinations of 
feddy are given with good accuracy by the Eddy Damped 
Quasi-Normal Markovian (EDQNM) two-point closure, 
previously known to reproduce well direct numerical sim- 
ulation results [13]. For helical flows, the EDQNM the- 
ory provides coupled equations for the energy and he- 
licity spectra [14], in which using (4) in an analogous 
way to [9] we find a very small correction of ^ e ddy that 
depends on the total amount of helicity and is of or- 
der Afeddy/^eddy ~ f3k max /a ~ 1CT 2 . Thus the pres- 
ence of helicity does not affect significantly the dissipa- 
tion at large scales and can be safely neglected in the 
eddy viscosity expressions. This eddy viscosity has a 
strong dependence in k and can also be obtained, in the 



limit /c/fc max — > 0, from the EDQNM eddy viscosity of 
Lesieur and Schcrtzer [15] using here an energy spectrum 
E(k) ~ fc 2 . The result reads 
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^oddy 



15A 



(7) 



with A = 0.36. The eddy viscosity f c ddy is thus an in- 
creasing function of time, see 2?th(*) in Fig. 2. 

The time-evolution of truncated Euler and Navier- 
Stokes spectra are compared in Fig. 3. At early times 
the value of is very small and therefore the NS vis- 
cosity vq is larger than ^ c ddy, as manifested by the NS 
dissipative zone in Fig. 3. a. As E t h(t) increases, both 
viscosities became equal (t = 2.7). Later, at t = 3.8, 
the Navier-Stokes spectrum crosses the truncated Euler 
one (Fig. 3b). The eddy viscosity f e ddy is then much 
larger than isq and the truncated Euler dissipative zone 
lies below the NS one, see Fig. 3c. This behavior is also 
conspicuous when the spectra are compared at maximum 
energy-dissipation time (t = 4.4 for truncated Euler and 
t = 5.6 for NS), see Fig. 3d. 
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FIG. 3: Compensated energy spectra of truncated Euler (• ■ ■ ) 
and Navier-Stokes ( x x x). a) t = 1.8, b) t = 3.8, c) t = 5.8. 
d) Maximum energy-dissipation time (t = 4.4 for truncated 
Euler and t = 5.6 for NS). 



The variation in time of f c ddy thus explains qualita- 
tively the different behavior of the truncated Euler and 
Navier-Stokes spectra in Fig. 3. We now proceed to 
check more quantitatively the validity of an effective dis- 
sipation description of thermalization in truncated Euler. 
To wit, we introduce an effective Navier-Stokes equation 
for which the dissipation is produced by an effective vis- 
cosity that depends on time and wavenumber. 
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We will use the effective viscosity obtained in [9] which 
is consistent with both direct Monte-Carlo calculations 
and EDQNM closure and is explicitly given by 



y eddyt 



-3.97fc/fc„ 



with v e ddy given in Eq. (7). 

We thus integrate Eq. (1) with the viscous term 
— v e ff(k)k 2 u a (h, t) added in the right hand side. The pa- 
rameter E t h that fixes the eddy viscosity in Eq. (7) is 
evolved using the effective NS dissipation by 



dE- 



th 



dt 



= £ 2v cS {k)k 2 E{k). 



(8) 



k=l 



This ensures consistency between the effective NS dissi- 
pated energy and the truncated Eulcr thermalized energy 
that drives f e ddy 

To initialize the effective NS equation we integrate the 
truncated Euler equation (1) with the initial condition 
(2) until the fc 2 -thermalized zone is clearly present (t = 
4.77). The value of E t h is then computed using equations 
(5). The low-passed velocity u<, defined by 



u<(r) = ^i(l + tanh[2(|k| 



fcth)])u k e t: 



is used as initial data for the effective Navier-Stokes dy- 
namics. 

Results of a truncated Euler and effective NS with 
fcmax = 85 are shown in Fig. 4. In Fig. 4. a the en- 
ergy and helicity dissipated in effective NS [E to t — E(t) 
and i/tot — H(t) respectively] are compared to i?th and 
-Hth showing a good agreement. Next, the temporal evo- 
lution of both energy spectra from the initial time t = 5.3 
(Fig. 4.b) to t = 20 (Fig. 4.c) is confronted, demonstrat- 
ing that the low-fc dynamics of truncated Euler is well 
reproduce by the effective Navier-Stokes equations. 

In summary, we observed the relaxation of the trun- 
cated Euler dynamics toward a Kraichnan helical abso- 
lute equilibrium. Transient mixed energy and helicity 
cascades were found to take place. Eddy viscosity was 
found to qualitatively explain the different behaviors of 
truncated Euler and (constant viscosity) Navier-Stokes. 
The large scale of Galerkin truncated Euler were shown 
to quantitatively follow an effective Navier-Stokes dy- 
namics based on a variable helicity-independent eddy vis- 
cosity. In conclusion, with its built-in eddy viscosity, the 
Galerkin truncated Euler equations appears as a minimal 
model of turbulence. 
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FIG. 4: Effective NS run with fc max = 85. a) Temporal evolu- 
tion of E t h ( — ), H t h ( — ) from truncated Euler, energy (. . .) 
and helicity (ooo) from effective NS. b-e) Temporal evolution 
of compensated energy spectra of truncated Euler (■ ■ ■ ) and 
effective Navier-Stokes (x x x) . 
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